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COMPENSATING POLARIZATION MODE DISPERSION IN FIBER OPTIC 

TRANSMISSION SYSTEMS 

RELATED APPLICATIONS 

This application is a continuation-in-part of U.S. Patent Application Serial No. 
5 09/505,313, filed February 16, 2000 and under 35 USC §1 19(e)(1), claims the benefit of U.S. 
Provisional Application Serial No. 60/243,522, filed October 23, 2000, which is incorporated 
herein by reference. 

BACKGROUND OF THE INVENTION 

This invention relates generally to fiber optics and more specifically to an apparatus 

1 0 for reducing the polarization mode dispersion in a fiber optic transmission system. 

Polarization mode dispersion (PMD) refers to variations in the time delay of a polarized 
optical signal traveling through an optical transmission system, such as a single-mode optical 
fiber. PMD arises in an optical fiber as a result of asymmetries in the optical fiber's core, 
such as core ellipticity created during optical manufacturing and bending stresses formed 

15 during installation. The asymmetries of the fiber's core cause random changes in the state of 
polarization (SOP) of optical signals propagating through the fiber. Different SOPs 
propagate through the optical fiber core at different relative speeds, e.g., some SOPs travel 
faster and some travel slower, resulting in a pulse width distortion of a transmitted optical 
signal relative to an input optical signal. Additionally, the asymmetries of the fiber's core are 

20 highly susceptible to environmental fluctuations, such as temperature or movement of the 
fiber, which occur as fast as milliseconds and result in a time varying pulse width distortion 
of the transmitted optical signal. 

The varying pulse width distortion is mathematically represented by a time delay 
between two orthogonally polarized principal states of polarization ("PSP") which form a 

25 convenient basis set to describe and characterize each SOP, and evaluate the effects of PMD 
in the fiber. Using the PSPs as a basis set, each SOP propagating through an optical fiber is 
represented by a linear combination of the two orthogonally polarized PSPs. The varying 
pulse width distortion of the SOP is a function of a varying delay between the PSPs. 
Theoretically, each PSP experiences a time of flight difference through the optical fiber, 

30 commonly known as differential group delay, resulting in a time delay between the two PSPs 
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at the fiber output. The output SOP is represented by a linear combination of the PSPs which 
are time delayed with respect to each other. A greater time delay between the PSPs 
corresponds to a larger relative difference between the input SOP pulse width and the output 
SOP pulse width. See for example CD. Poole and R.E. Wagner, "A Phenomenological 

5 Approach to Polarization Dispersion in Long Single-Mode Fibers." Electronic Letters, Vol. 
22, pp. 1029-1030, September 1986, which is incorporated by reference herein. 

Optical fibers have a differential group delay (DGD) between the two PSPs on the 
order of 0.1 ps/km. In older fiber optic cables, such as the cables used in terrestrial networks, 
the DGD is on the order of 2.0 ps/km and results in time delays of about 50 picoseconds for 

1 o transmission distances of only several hundred kilometers. As the demand for faster optical 
data transmission increases, such as from gigabits per second to terabits per second, optical 
pulse width distortion due to PMD will become one of the factors limiting data transmission 
rate. 

SUMMARY OF THE INVENTION 

15 A real-time optical compensating apparatus reduces first-order PMD in an optical 

fiber by determining the PSPs of the optical fiber and delaying one PSP with respect to the 
other. 

In one aspect, the invention features an optical compensating apparatus including a 
polarization controller configured to receive an optical signal propagating through an optical 

20 medium, determine the principal states of polarization of the optical medium, and transform 
the optical properties of the optical signal based on the principal states of polarization. 

Embodiments of this aspect may include one or more of the following. The 
polarization controller includes a polarimeter configured to receive and to measure the state 
of polarization of the optical signal, and a polarization transformer arranged in the path of the 

25 optical signal after the polarimeter and configured to transform the polarization of the optical 
signal. The polarization controller includes a central processing unit programmed to monitor 
the time averaged state of polarization of the optical signal and determine the principle states 
of polarization of the optical medium based on the time averaged state of polarization. The 
central processing unit analyzes time averaged states of polarization via the relationship, 
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DOP = |^|= ' — ' , to determine the principle states of polarization by calculating 

\r k -(a-r k )a\ 

the vector a . The central processing unit is programmed to determine a magnitude of time 
delay between the principle states of polarization based on the time averaged state of 
polarization. The central processing unit is programmed to determine the relative order, in 

5 time, of the principle states of polarization. The compensating apparatus includes a delay 
controller arranged in an optical path of the optical signal after the polarization controller. 
The delay controller is configured to receive the transformed optical signal and to generate a 
signal proportional to the PMD time delay between the principle states of polarization. The 
delay controller includes a DOP sensor, such as a polarimeter or spectrum analyzer. The 

10 delay controller includes a delay transformer arranged in the path of the transformed optical 
signal before the polarimeter. The compensating apparatus further includes a PMD 
generating module arranged in the optical path of the medium before the polarization 
controller and configured to receive the optical signal from the optical medium and generate 
a known PMD. The compensating apparatus further includes a calibrating unit arranged in 

15 an optical path of the optical signal after the polarization controller. The calibrating unit is 
configured to receive the transformed optical signal and determine the transformation of the 
polarization controller. The calibrating unit includes at least one detector for detecting one or 
more components of the optical signal. 

In another aspect, the invention features, an optical compensating apparatus for 

20 reducing PMD in an optical signal transmitted through an optical medium. The optical 
compensating apparatus includes a polarization controller configured to receive an optical 
signal propagating through the optical medium, determine the principal states of polarization 
of the optical medium and a magnitude of time delay between the principal states of 
polarization; and a delay controller arranged in an optical path of the medium after the 

25 polarization controller and configured to receive the transformed optical signal, wherein the 
polarization controller transforms the optical signal based on the principal states of 
polarization and the delay controller compensates PMD of the optical signal by reducing the 
time delay between the principle states of polarization based on the magnitude of time delay. 
Embodiments of this aspect may include one or more of the following. The 

30 polarization controller includes a polarimeter. The polarization controller includes a 
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polarization transformer arranged in the path of the optical signal after the polarimeter. The 
delay controller is configured to generate a signal proportional to the PMD time delay 
between the principle states of polarization. The delay controller includes a DOP sensor, 
such as a polarimeter or spectrum analyzer. The delay controller includes a delay 

5 transformer arranged in the path of the transformed optical signal before the polarimeter of 
the delay controller. The delay controller includes a delay transformer arranged in the path 
of the transformed optical signal after the polarimeter of the polarization controller and 
before a polarimeter of the delay controller. 

In another aspect, the invention features an optical compensating apparatus for 

10 reducing PMD in an optical signal transmitted through an optical medium. The apparatus 
includes a polarization module configured to receive an optical signal propagating through 
the optical medium, determine the principal states of polarization of the optical medium, and 
generate a signal for transforming the polarization of the optical signal; and an optical 
transformer arranged in an optical path of the medium after the polarization module and 

15 configured to transform the optical signal and reduce a time delay between the principle 
states of polarization based on the signal received from the polarization module. 

Embodiments of this aspect may include one or more of the following. The 
polarization module includes a polarimeter. The optical transformer includes a polarization 
transformer and a delay controller. The polarization transformer and a delay controller each 

20 include a polarimeter. The polarimeter of the delay controller is configured to generate a 

signal proportional to the PMD time delay. The delay controller includes a delay transformer 
arranged in the path of the transformed optical signal before the polarimeter of the delay 
controller and after the polarization transformer. 

In another aspect, the invention features a method of reducing PMD of an optical 

25 signal propagating in an optical medium. The method includes determining the principal 
state of polarizations of the optical medium with a polarization controller, and transforming 
the polarization of the optical signal with a polarization transforming device based on the 
polarization of the first principal state of polarization. 

Embodiments of this aspect may include one or more of the following. The method 

30 further includes determining the time delay between a first principal state of polarization and 
a second principal state of polarization. The method further includes delaying the first 
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principal state of polarization with respect to a second principal state of polarization. The 
method further includes delaying a first principal state of polarization with respect to a 
second principal state of polarization. The method further includes determining the relative 
order, in time, of the principle states of polarization. The method further includes monitoring 
5 a time averaged state of polarization of an optical signal propagating through the optical 
medium. The step of determining the principle states of polarization includes analyzing the 
time averaged state of polarization. The time averaged states of polarization follow the 

relationship, DOP = \f k \ = LJJL! ? an( i determining the principle states of 

— > 

polarization includes calculating the vector a . The vector a is determined by analyzing r k 
10 with nonlinear curve fitting techniques or matrix methods. The method further includes 
calibrating a polarization controller for transforming the optical signal to determine a 
plurality of control settings for the polarization controller by monitoring which control 
signals when applied to the polarization control transform an input optical signal having a 
specific state of polarization into an output signal having different state of polarization. 
1 5 In another aspect, the invention features a method of compensating PMD of an optical 

signal propagating through an optical medium. The method includes recording the time 
averaged state of polarization of an optical signal propagating through an optical medium; 
and analyzing the time averaged state of polarization of the optical signal via a pertubative 

expansion of the time averaged state of polarization, r , to extract a plurality of parameters 
20 from the time averaged state of polarization that characterize the PMD of the optical signal, 
wherein the pertubative expansion is expressed as 

DOP 2 S |r \ 2 =\jSo/ +( l -^uise4M D )\SoJ + °{ Ao > 3 }- The plurality of parameters, for 

example, can be used to characterize the location of the principle states of polarization on the 
Poincare sphere, the magnitude of time delay between the principle states of polarization, and 
25 the frequency dependence of the PMD. 

The invention has various advantages including, but not limited to, one or more of the 
following. The apparatus for compensating PMD operates in real time and does not require a 
fast detector. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram of a PMD compensating apparatus; 

Fig. 2 is a block diagram of a polarimeter shown in Fig. 1 ; 

Fig. 3 is a block diagram of a delay module shown in Fig. 1 ; 

Fig. 4 is a schematic representation of a Poincare sphere; 

Fig. 5A is a schematic representation of an optical signal propagating through an 
optical fiber of Fig. 1; 

Fig. 5B is a graphical representation of the optical signal propagating through an 
optical fiber of Fig. 1; 

Fig. 6 is a schematic flow chart of the PSP algorithm; 

Fig. 7A is a schematic representation of the principal states of polarization of an 
optical signal propagating through an optical fiber of Fig. 1; 

Fig. 7B is a graphical representation of the degree of polarization as a function of the 
time delay, t, between and the relative power, a, of each principal state of polarization; 

Fig. 8 is a cross-sectional view of the Poincare sphere of Fig. 4 taken along the 
equator; 

Fig. 9 is a block diagram of another PMD compensating apparatus; 

Fig. 10 is a block diagram of another PMD compensating apparatus; 

Fig. 1 1 is a block diagram of a apparatus for calibrating a polarization transformer; 

and 

Fig. 12 is block diagram of an expanded view of a delay module integrated with a 
calibration unit. 

DESCRIPTION OF THE EMBODIMENTS 

Referring to Fig. 1, compensating apparatus 10 includes a polarization controller 100 
and a delay controller 200. Compensating apparatus 10, when placed between an output 21 
of an optical fiber 22 and an input 235 of optical receiver 240, reduces the PMD of optical 
signals transmitted by optical transmitter 15 though optical fiber 22. 

Polarization controller 100 includes a lens 104, a beam splitter 105, a polarimeter 
110, and a polarization transformer 108. Lens 104 positioned at an input 102 of polarization 
controller 100 collimates optical signals (not shown) from output 21 of optical fiber 22 along 
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an optical path 160. Optical path 160 extends from input end 102, through beam splitter 
105, polarization controller 108, and out output end 103. Beam splitter 105 redirects a 
portion of the optical signal propagating along beam path 160 into polarimeter 110 which 
detects the redirected optical signals and sends a series of electronic signals through cables 

5 122 to an I/O port 1 19 of a computer 120. Computer 120 uses the electronic signals in an 
algorithm stored in the computer's CPU to determine the principal states of polarization 
(PSPs) of optical fiber 22 and sends control signals to modify the settings of a first retarder 
140 and a second retarder 150 in polarization transformer 108. First retarder 140, e.g., a 
quarter-waveplate, and second retarder 150, e.g., a half-waveplate, transform the polarization 

10 of the PSP such that light exiting polarization controller 100 is linearly polarized and aligned 
to the x- and y- axis of delay controller 200. 

Delay controller 200 includes a delay module 170, a beam splitter 165, a polarimeter 
210, a controller 220, and a mirror 202. An optical beam path 162 extends between an input 
161, through delay module 170, and beam splitter 165. Mirror 202 reflects optical signals 

15 out of delay controller 200 through output 163 and into input 235 of receiver 240. Beam path 
1 62 at input 161 is collinear with beam path 160 from polarization controller 100 such that 
collimated optical signals exiting output end 103 propagate along beam path 162. 

After the optical signals pass through delay module 170, beam splitter 165 redirects a 
portion of the optical signal propagating along beam path 162 into polarimeter 210. 

20 Polarimeter 210 detects the redirected optical signals and sends a series of electronic signals 
via cables 215 to a control circuit 220. Polarimeter 210 is similar in structure to polarimeter 
110 described below. Control circuit 220 uses the electronic signals sent from polarimeter 
210 to determine the time delay between the PSPs in optical fiber 22 and then sends a control 
signal via cable 130 to delay controller 170. The control signal modifies the settings of delay 

25 controller 1 70 such that the time delay is reduced between the two PSPs transmitted through 
outlet 21 of optical fiber 22. 

Referring to Fig. 2, polarimeter 110 includes three beam splitters 1 14, 1 16, 1 17, and a 
mirror 1 19 spaced along an optical beam path 112. Beam splitters 1 14, 1 16, 1 17, and mirror 
119 couple optical signals propagating along beam path 1 12 towards detector modules 1 14a, 

30 1 16a, 1 17a, 1 19a, respectively. Detector module 1 14a includes a detector 1 14b for 

measuring the total power of an optical signal. Detector module 1 16a includes a polarizing 
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beam splitter 1 16b and a detector assembly 1 16c having a first detector 1 16d and a second 
detector 1 16e. Similarly, detector module 1 17a includes a polarizing beam splitter 1 17b and 
a detector assembly 1 17c. Detector module 1 19a includes a polarizer 1 19f, e.g., a quarter- 
waveplate, a polarizing beam splitter 1 19b, and a detector assembly 1 19c. Each detector 

5 module measures specific optical properties of the optical signal and sends an electronic 
signal proportional to each measured property to computer 120 via cables 122. 

Referring to Fig. 3, delay module 170 includes an input polarizing beam splitter 171, 
an optical delay assembly 174, and an output polarizer 172. Polarizing beam splitter 171 
separates the two PSPs transmitted through optical fiber 22 and polarization controller 100 

1 o such that PSP 1 , delayed with respect to PSP2, propagates along a fixed optical path 1 73 to 
polarizing beam splitter 172, and PSP2 propagates along a variable optical path 175. 
Variable optical path 175 includes optical delay assembly 174, e.g., a translatable (Arrows) 
mirror, which delays PSP2 with respect to PSP1. PSP1 and PSP2 recombine in polarizing 
beam splitter 172 and continue propagating along beam path 162. 

1 5 Referring to Fig. 4, a convenient and intuitive graphical representation of SOPs is a 

Poincare sphere 500. A SOP is defined in terms of a Stokes vector of Formula 1 : 

{S 0 S, S 2 S 3 ) (1) 

where 

20 So = E 2 x + E 2 y (2) 

S!=E 2 X -E 2 y (3) 

5 2 = 2E x E y cos(0) (4) 

5 3 = 2E x Eysin(0) (5) 

and E x and E y are the magnitudes of the x and y component electric field complex 
25 amplitudes, respectively. 0 is the relative phase between the two. The parameters si, s 2 and 
s 3 are represented by the relationship s\= (Si)/(S 0 ), where i is 1, 2, or 3, and can be used to 
convert the Stoker parameters to corresponding x, y, and z components in a three 
dimensional Cartesian coordinate system. In a three dimensional Cartesian coordinate 
system, Poincare sphere 500 is defined by a set of points containing all possible SOPs. As 
30 seen in Fig. 4, all linear polarization states are located on an equator 5 10 of sphere 500, while 
left and right circular polarizations are located at a north pole 520 and a south pole 530, 
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respectively. All other points represent elliptical polarizations which lie away from equator 
5 1 0 and poles 520, 530. Each SOP on sphere 500 is identifiable by its latitude 2w and 
longitude 2X by using equations: 

si = cos(2w)cos(2A,) (6) 
5 s 2 = cos(2w)sin(2X) (7) 

S3 = sin(2w) (8) 

where any two orthogonal SOPs lie directly opposite each other, e.g., a linear vertical 
polarization 512 is 180 degrees away from a linear horizontal polarization 514 on equator 
510. 

10 In operation, transmitter 15 sends a polarized optical signal to polarization modulator 

20 which modulates the state of polarization ("SOP") of the optical signal, e.g., from vertical 
to right circular to elliptical, with a frequency of about 10 kHz to about 100 MHz. The 
modulating frequency is fast enough to measure and compensate varying PMD on a 
millisecond timescale. The modulating frequency is limited by the response time of the 

1 5 detectors used in the polarimeters. 

Referring to Fig. 5 A, as an input SOP 550 propagates from polarization modulator 
(not shown) through optical fiber 22 towards output 21, the SOP of the signal randomly 
changes. Each SOP propagates through the optical fiber at a different speed, e.g., some SOPs 
travel faster and some travel slower, resulting in a varying pulse width distortion of an optical 

20 signal 560 at output 21 . To a first order approximation, optical fiber 22 has two discrete 
group delays, one for each of two orthogonal PSPs, i.e., PSP 562 and PSP 564. Referring to 
Fig. 5B, optical signal 560, a SOP, is a linear combination of PSP 562 and PSP 564. A time 
delay 565, e.g., 40 ps, between PSP 562 and PSP 564 creates PMD in the output signal, i.e., 
the width of signal 560 is greater than signal 550. 

25 At outlet 21, optical signal 560 propagates into compensating apparatus 10 which 

reduces the PMD in signal 560 with polarization controller 100 and delay controller 200. 
Polarization controller determines the PSPs of the optical fiber 22 and transforms the two 
PSPs to x and y linearly polarized states aligned with the x and y optical axis of delay module 
170. Delay controller 200 measures the time difference between the two transformed PSPs 

30 and delays one PSP relative to the other. 
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After exiting optical fiber 22, optical signal 560 travels through polarization 
controller 100 along beam path 160, until beam splitter 105 redirects a portion, e.g., about 
1%, of the optical signal into polarimeter 1 10 for analysis. The amount of optical signal 
redirected into polarimeter is sufficient such that the redirected optical signal is measurable 

5 by each of the detectors in polarimeter 1 10. 

Referring back to Fig. 2, beam splitters 1 14, 1 16, 1 17 evenly divide the optical signal 
entering polarimeter 110 into four separate optical signals. Detector module 1 14a measures 
the power of the first optical signal, i.e., So. Detector module 1 16a analyses the second 
signal by measuring the difference between optical signals having polarization components 

10 oriented in the x and y direction, i.e., Si = E 2 X - E 2 y . Optical signals having only x oriented 
polarization components result in a measurement of +1 by detector module 1 16a, and optical 
signals having only y oriented polarization components result in a measurement of -1. 
Detector module 1 17a analyzes the third signal by measuring the difference between optical 
signals having polarization components oriented 45 degrees with respect to the x and y 

15 direction, i.e., S2 = E 2 +45 = 2ExE y cos(0) (where 0 is the phase between E x and E y ). Detector 
module 1 1 9a analyzes the fourth signal by measuring the difference between optical signals 
having left and right circular polarization components, i.e., S3 - 2E x E y sin (0) where 0 is the 
phase between E x and E y . Optical signals having only right circular polarization results in a 
measurements of +1 , and left circular polarization results in a measurement of -1 . 

20 Polarimeter 110 measures each of the Stokes parameters, Equations 2-5, and sends electronic 
signals proportional to each measurement to computer 120. 

Referring to Fig. 6, computer 120 receives the electrical signals from polarimeter 110 
and runs an algorithm 600 stored in the computer's CPU to determine the location of the SOP 
on a Poincare sphere (S10), to calculate a degree of polarization (DOP) for each SOP (S20), 

25 to determine the two orthogonal PSPs for the optical fiber (S30), and to calculate a 

polarization transformation which converts the PSPs of the fiber to linear x and y polarization 
states aligned with the x and y optical axis of delay module 170 (S40). 

The CPU determines the location of the SOP on the Poincare sphere by relating each 
of the electrical signals from the polarimeter to its corresponding Stokes parameter, equations 

30 6-8, and then calculating w and A,, i.e., the coordinates of the SOP on the Poincare sphere. 
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Each of the Stokes parameters is also used by the CPU to calculate degree of polarization 
(DOP). The ratio of Stokes parameters shown below 

DOP 1 = S ' + S \ + S ' (9) 

is less than or equal to unity. 
5 Referring to Fig. 7A, two PSPs, i.e., square shaped pulses of length To, have a first 

order PMD delay of x and a relative power with respect to each other described by the 
parameter a. To a first order approximation, i.e., when t/tq is less than 1, the DOP of the 
two square PSPs, equation 9, can be rewritten as a function of delay, x, and relative power, a: 

DOP(x,a) = [1 + 4 a (x/x 0 ) ((x/x 0 )-2)(l- a)] 1/2 . (10) 

1 o A plot of Equation 1 0 (Fig. 7B) graphically shows how DOP depends on delay and relative 
power of the PSPs. At constant delay, x, DOP is at a minimum when both PSPs have equal 
power, whereas DOP is at a maximum when only one PSP has all of the power, i.e., a is 1 or 
0, respectively. At constant relative power, a, DOP is inversely related to the delay between 
the two PSPs. At a constant time delay between the two PSPs, the DOP depends on the SOP. 

15 When the SOP is a 50/50 mixture of both PSPs, i.e., each PSP has equal power, the DOP will 
be at a minimum, whereas the DOP will be unity when the SOP is aligned with one PSP, i.e., 
one PSP has all the power. As the value of x/x 0 approaches 1, the first order approximation 
of PMD fails and equation 10 is no longer valid. 

Referring to Fig. 8, the linear polarization states represented by circle 570 are 

20 synonymous with equator 510 of the Poincare sphere 500 (see Fig. 4). Assuming that optical 

fiber 22 includes an x-horizontal linear PSP 575 and the y-vertical linear PSP 577, i.e., two 

orthogonal PSPs, all other points on the circumference of circle 570 represent linear states at 

different orientations. A SOP 578 represents one possible linear SOP of an optical signal of 

optical fiber 22. SOP 578 contains components of both x and y polarizations, i.e., SOP 578 

25 is a weighted linear combination of PSP 575 and PSP 577. Depending on the amount of time 

delay between PSP 575 and PSP 577, SOP 578 has a DOP that is less than or equal to unity. 

As the angular distance, 2k, of SOP 578 to a PSP decreases, DOP increases. At the critical 

angular distance, 2k cnh a SOP 579 is equally distant from PSP 575 and PSP 577, i.e., SOP 
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10 



15 



20 



579 is a 50/50 mixture of PSP 575 and PSP 577, resulting in a minimum DOP. The 
definitions of the Stokes parameters, such as Si, provide a relation between a and the radial 
distance 2X such that Equation 10 takes the form. 



DOP(x,2X) = 



1 + 



-2 



1 + 



o V T o 



COS 



(21) 



(11) 



In general, the DOP is a function of 2X\ the angular distance between the SOP and any PSP 
on the Poincare sphere. 2V is a function of the longitude distance, 2X, and the latitude 
distance, 2w. Using the definitions of the Stokes parameters, Equation 10 is rewritten in the 
form. 



DOP(x,2V) = 



1 + 



-2 



1+ 



--2 



COS 



(2V) 



(12) 



If the PSPs of an optical fiber are unknown, computer 120 runs algorithm 600 while 
polarization modulator 20 changes the SOP. For each SOP, algorithm 600 calculates both 
the location of the SOP on the Poincare sphere and the DOP. The algorithm sends these 
values into memory and repeats the cycle. Algorithm 600 stops collecting data points and 
fits, e.g., by linear-least-squares, the data in memory to find the maximum DOP, i.e., a DOP 
of unity corresponds to a SOP which represents one of the PSPs. Algorithm 600 collects 
enough data points so that the DOP as a function of SOP is well represented. Collecting too 
few data points leads to incorrect fitting results, whereas collecting too many data points is 
time consuming and allows environmental changes, i.e., temperature and stress on the fiber, 
to affect the location of the PSPs on the Poincare sphere. Once the CPU calculates the 
identity of the PSPs, algorithm 600 calculates a polarization transformation necessary to 
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transform the PSPs into linearly polarized x-horizontal and y-vertical PSPs aligned to the x- 
horizontal and y-vertical axis of delay module 170. 

Algorithm 600 uses Stokes parameters and Jones matrices representing the PSP 
polarization states, retarder 140, and retarder 150 to calculate the settings of the retarders 
5 which transform the PSPs into linear horizontal and vertical states. Algorithm 600 begins 
with an arbitrary polarization state A 



8 



(13) 



10 



15 



20 



25 



where S 2 + e 2 = 1, and a horizontal polarization state, i.e., a x-horizontal state, is represented 
by the Jones matrix, 

T 



x - horizontal - 



0 



(14) 



The angle, A, between a major axis of the arbitrary state and the axis of a Cartesian 
coordinate system is 

A = 2 tan 4 [(2 8 s cos(<p))/8 2 - s 2 )]. (15) 

Algorithm 600 converts the known PSPs from Cartesian coordinates into Jones matrices, and 
determines A. 

Once A is known, Algorithm 600 calculates how to orient retarder 140 to convert the 
PSPs into linearly polarized states having an angle, A, between the x-axis of the Cartesian 
coordinate system and the axis of the linearly polarized states. Next, Algorithm 600 
calculates how to orient retarder 150 to rotate, i.e., by p, the linearly polarized states so that 
they coincide with x-horizontal and y-vertical polarization states. 

A complete transformation of an arbitrary state into a linearly x-horizontal 
polarization state using a quarter- waveplate as retarder 140 and a half-waveplate as retarder 
150 is 



[R(-p/2) X HW X R(p/2)] X [R(-A) X QW X R(A)] 



(16) 



where 
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R(D) = 



cos(D) sin(D) 
-sin(£>) cos(Z))J' 



(17) 



HW = 



i 0 
0 -i ' 



and 



(18) 



5 



QW = 




0 



(19) 



0 e 



The computer sends electrical control signals, based on the transformation 
calculations, both to retarder 140 to transform the PSPs to linear PSPs and to retarder 150 to 

1 o rotate the linear PSPs to the x-horizontal and y-horizontal axis of the delay module. 

The transformed optical signal propagates into delay module 170. Referring back to 
Fig. 3, polarizing beam splitter 171 redirects x-horizontally polarized optical signal, i.e., PSP 
2, through a variable delay line and y- vertically polarized optical signal, i.e., PSP 1, through 
a fixed delay line. Polarizing beam splitter 172 recombines the x-horizontally and y- 

1 5 vertically polarized optical signals. Before the recombined signal exits delay controller 200, 
beam splitter 165 redirects a portion of the optical signal into polarimeter 210. Polarimeter 
210 is similar to polarimeter 1 10 described above and sends electrical signal proportional to 
the stokes parameters to control circuit 220, 



20 signal As shown in Equation 10, DOP is a function of both the time delay, t, and the SOP of 
the optical signal, i.e., the relative power, a, of each PSP. Control circuit 220 calculates 
DOP and determines a time average DOP as polarization modulator 20 modulates the SOP of 
the input optical signal and the relative power, a, of each PSP. Referring back to Fig. 7b, the 
time averaged DOP, i.e., as a changes between values of 0 and 1, is at maximum for zero 

25 delay between PSP1 and PSP2, whereas the time averaged DOP decreases as the delay 

between the PSPs increases. Control circuit 220 sends electrical signals to delay assembly 



The control circuit, e.g., a microprocessor, calculates the DOP of the recombined 
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1 74 to adjust the time delay between PSP1 and PSP2 such that the time averaged DOP is 
maximized. 

Control circuit 220 averages the DOP for a time period that is sufficient to 
characterize the DOP for several SOPs. The shortest DOP averaging time period is set by the 

5 rate at which polarization controller 100 aligns the PSPs to the x- and y-axis of the delay 
module. At longer time periods, collecting too many data points is time consuming and 
allows environmental changes, i.e., temperature and stress on the fiber, to affect the location 
of the PSPs on the Poincare sphere. 

In an alternative embodiment, referring to Fig. 9, a CPU 720 of a compensating 

10 apparatus 700 receives electrical signals from a polarimeter 730 and a control unit and runs 
an algorithm, as will be described below, to determine the location of the two PSPs of optical 
fiber 740 on the Poincare sphere, and the magnitude of PMD time delay between and the 
relative order, in time, of the two PSPs. The Control unit includes a DOP sensor 770, such as 
a polarimeter or a spectrum analyzer, and a CPU 780. DOP sensor 770 monitors a portion of 

15 the optical signal exiting a delay module 760 and sends electric signals to CPU 780. CPU 
780, in turn, calculates a value of DOP, DOPsign, based on the electric signals sent from 
DOP sensor 770 and sends an electrical signal proportional to DOPsign back to CPU 720. 
CPU 720, in turn, uses this value of DOPsign to determine the relative order, in time, of the 
two orthogonal PSPs, e.g., which PSP leads and which PSP trails. 

20 During operation, CPU 720 determines which PSP is leading and which is trailing by 

comparing two different DOPsign values, one DOPsign value determined when the delay 
module retards one PSP and the other DOPsign value determined when the delay module 
retards the other PSP. A higher DOPsign value corresponds to less PMD time delay and a 
correctly retarded PSP. When the wrong PSP is delayed, i.e., the trailing PSP is further 

25 delayed in time relative to the leading PSP, the PMD time delay increases and the DOPsign 
value decreases. Typically, CPU 720 records and compares two values of DOPsign at a rate 
faster than the rate at which PMD changes in the optical fiber. After determining which PSP 
is leading and which is trailing, CPU 720 continuously monitors the DOPsign value 
calculated by CPU 780 to determine whether or not the leading and trailing PSPs have 

30 switched order. If the PSP order switches, the DOPsign value will decrease because the 
delay module will further delay the "new" trailing PSP thereby generating more PMD time 
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delay between the two PSPs. When CPU 720 detects a decrease in the DOPsign value, the 
CPU uses this information to generate electric signals which when sent to a polarization 
transformer 750 cause the transformer to rotate the PSPs relative to the delay module such 
that the delay module will retard the "new" trailing PSP. Alternatively, without rotating the 
5 polarization of the PSPs, CPU 720 can send an electric signal to the delay module causing 
the module to change which PSP is retarded, such as by changing the optical path delay. See 
Fig. 3. 

CPU 720 runs an algorithm that records electric signals from the polarimeter and the 
control unit and stores these values into memory. The algorithm processes the stored values 

1 0 from the polarimeter to determine where the PSPs of the optical fiber are located on the 

Poincare sphere and the magnitude of the time separating the two PSPs. The algorithm also 
processes the stored values from the control unit, as discussed above, to determine the 
relative order of the PSPs. Once CPU 720 has determined the location of the PSPs and 
identified the leading and trailing PSP, CPU 720 calculates a polarization transformation that 

15 will orient the PSPs relative to a delay module 760 such that the module will retard the 
leading PSP relative to the trailing PSP. CPU 720, in turn, sends electrical signals to 
polarization transformer 750 to produce the calculated polarization transformation. CPU 720 
determines which electric signals will generate the desired transformation by using a lookup 
table containing input polarizations, values of applied electric signal voltages, and output 

20 polarizations. The lookup table can be preprogrammed by a system user into CPU 720 prior 

to operation or determined via a calibration procedure described below. After determining 

the magnitude of the PMD time delay, CPU 720 sends electrical signals to delay module 760 

to set the overall magnitude of time delay between the two PSPs. 

CPU 270 processes each of the signals from the polarimeter to generate a time 

—> 

25 averaged state of polarization of a stokes vector, r . CPU 270 uses the time averaged state of 

polarization to extract an estimate of the location of the PSPs and the magnitude of their 

— > 

separation via a pertubative expansion of r . 

DOP 2 m\r\ 2 =\$ 0I \ 2 Hl-Aco 2 puL y mD )\s 0 J+o{^ 3 } ( 20 > 

where r is the time averaged state of polarization of a stokes vector measured by the 

30 polarimeter, s oj| and s 0,3 are two components of the SOP, o pu i se is the signal bandwidth, 
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t 2 pmd is the magnitude of the PMD time delay and can be expressed as |Q| 2 , and 0{Ag) 3 } 
represents higher-order terms of the pertubative expansion. Equation 20, minus the higher- 
order terms, 0{Aoo 3 }, is equivalent to equation 10 and is generated by starting with a two- 
mode, pure-polarization input field, written aft), aft) is transformed by an optical fiber into 
an output field bft), expressed as 

b{t) = T(w)a(t) = ^f( w )e^T(w)a 0 (21) 

where /is the signal spectrum and a 0 is a complex 2-dimensional vector. The time averaged 
state of polarization is represented by 



f = Jdtb'(t)ab(t) 



(22) 



where a is a vector of Pauli matrices. Using equations 21 and 22, the time averaged state of 
polarization can be expressed as 

7 ■ ' f^W^ 5 S iff 1 f ( ^0 )| 2? ( <^, ) (23) 

PMD is characterized by a (Stokes-space) PMD vector, Q , which can be defined by the 
relation 

do w w (24) 

By using equation 24, and by defining a center of the signal spectrum coo as 

J^I/(^)| 2 (^"^o)-0 
J in 

and a nominal SOP as 

s 0 sa;r'(ffl 0 )ff r(^ 0 ) = o, 

a perturbative expansion for time averaged state of polarization can be expressed as 

Q^x^o + Q x(Q x^ 0 ) J+0(A&> 3 } (25) 



-» Acq 2 , 

t - J * w pulse 

r ™ so- 



2 

where the signal bandwidth is characterized by Aco 2 pulsc = ^dw\ J(w)| 2 |(w - w 0 ) 2 / 2;r . 
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The Qxflx term can be described physically as a reduction of components orthogonal to Q 
due to first-order PMD and the £2 w x so term can be described physically as a rotation due to 
second-order PMD. Rearranging equation 25 yields 

1 A 2 2^,1 

— A© nilUe T ^ s 0 ,±-f - 

2 P ulse PMD u 2 

Using equation 26, the physical effects of PMD of the optical fiber on an SOP can be 
described physically as generating a rotation and a rescaling of the SOP: 



f = s 0 -^AoV^t ^s 0 ,±+^A©^Q ffl x s 0 +o{Aco 3 (26) 



r = 



Rotate (Qj] [scale (Q)] s Q . (27) 



In terms of first-order PMD, a rescaling means that DOP is a function of the angle 

-> 

between s o and Q . Although discussed in terms of first- and second- order PMD, the 
10 perturbative expansion for time averaged state of polarization can be used to develop 
analytical models for higher-order PMD which can be used to extract higher-order PMD 
parameters from DOP measurements. For instance, a CPU in a PMD compensating 
apparatus capable of compensating higher-order PMD can use the higher-order PMD 
parameters derived from equation 20 to generate control signals to reduce higher-order PMD 
1 5 in the optical signal. 

For simplicity, we limit our discussion to first-order PMD calculations. Given a 
number, N, of time-averaged SOP measurements r 1} ...,r k , resulting from unknown and 
varying input polarization input fields ao, the samples, n ideally obey the relation 

DOP = \f k \= (28) 
\r k -(a-r k )a\ 

20 for first-order PMD, where a is an unknown vector to be estimated by the CPU. In Stokes 



space, 



rk 



forms the surface of an ellipsoid within the Poincare sphere whose major axis is 



aligned with the PSP directions +a and - a. Vector a satisfies the relationship (1-| a | 2 ) 2 = 1 - 

— > 

A© puise x P md which relates the magnitude of a to the PMD time delay between the PSPs. 
The CPU can estimate the vector a by mathematical methods such as eigenvalue 
25 calculations or nonlinear curve-fitting. In nonlinear curve-fitting, the CPU records a series of 



DOP data, and coverts them into magnitudes, <& □ 



rk 



and directions f k = r/ 



. The 
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CPU performs a nonlinear fit by minimizing the distance between the measured DOP (dk) 
and the expected DOP (f(r k, a )). For example, the least-squares estimate of a is 

^.ri = argmini^(&) 
b 

where 

5 W>lK-/(^)] 2 (30) 

k 

f(r k b) = H , b|2 - (31) 
' l? k -(b-r k )b| 

where the function f( r k, b ) gives the theoretical degree of polarization (DOP), or 

-» 

magnitude of r . 

— > 

In eigenvalue calculations, the CPU can estimate a by using matrix method 
1 0 techniques to maximize the second moment of the points on an unknown projection. 

Typically, CPU rescales the measured DOPs, r kj before performing moments calculations 
to artificially sharpen the contrast in DOP. CPU rescales the time-averaged SOP 
measurements via the relationship 

q k =S(|f k |)f k (32) 
1 5 where S(|r k |) is equal to ((DOP k -DOP min )/ 1 -DOP mi „). CPU uses the scaled vectors, q k , to 
calculate an estimate of the PSP axis 

««*2 =argminF 2 (£) 
\b\=\ 

where 

f 2 (*) = £(H) 2 < 34 > 

k 

20 The function F2 generally reaches a maximum near the direction of the PSPs, because 

PMD tends to pull the SOP towards this axis. CPU can also use the unsealed time-averaged 

DOP measurements, r , to estimate a , but large values of DOP m i n decrease the contrast 
between the measurements which, in turn, can decrease the precision of the method because 
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random nonuniformities in the input polarizations can dominate the function F 2 and reduce 

— > 

the accuracy of determining a . Equation 34 can be rearranged to provide 

F 2 (b) = b T b^b T Mb (35) 

which reduces to the very simple form: 
5 Ma est2 ^ka est2 (36) 

equation 36 can be solved using Lagrange multipliers. The estimated a is the eigenvector 
corresponding to the largest eigenvalue of the 3x3 matrix M. The magnitude of the PMD 

time delay is equivalent to the magnitude of the vector a . The CPU determines the 

magnitude of a by analyzing the statistics, such as mean and minimum, of the DOP data. 
10 Random fluctuations in the density of sampled SOP's can effect the precision at 

— > 

which CPU determines the estimated a . To increase precision, the CPU can record more 
samples. Recording more samples reduces the overall speed at which the compensating 
apparatus can reduce PMD in the optical signal. Alternatively, the CPU can increase 
precision by correcting random fluctuations in the density of a smaller subset of recorded 
1 5 samples. CPU reduces the random fluctuation by using a moment function to remove DOP: 

D(b) = b T Wr/Ab ^ b T M 2 b (37) 

\k J ) 

The matrix method also has an advantage that much of the calculation can be 
performed in real-time as the CPU records data from the polarimeter. Each data point can be 
used to increment a matrix, M, with very few floating operations by continuously updating M 
20 to include new data and throw out old data. The CPU can update Musing the relationship: 
M ( * ) =(l»tf)Af ( *" ,) +^ 9 ; (38) 
When using this scheme, the CPU can maintain an updated Mat all times so that any 
calculations, such as an eigenvalue decomposition to estimate the PSPs, can keep pace with 
data acquisition. 
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The vector, a , estimated by the CPU has three degrees of freedom, one degree for a 

and two degrees for the direction of a . The relationship between the 

-> -» 
magnitude of vector a and the PMD vector, Q , which also has three degrees of freedom, is 

Q = ±(2|a| 2 /< ls J' 5 a. (39) 

5 where a) 2 pu i se is the bandwidth of the optical signal and can be determined by measuring the 
pulse shape of the optical signal with a spectrum analyzer. The CPU can be programmed 
with an optical signal bandwidth that was measured prior to operation because changes in 
PMD generally do not affect bandwidth. Alternatively the compensating apparatus can 
include a spectrum analyzer which monitors the bandwidth of the optical signal and sends 

1 0 electronic signals proportional to that bandwidth to the CPU. 

When calculating an estimated PMD vector with equation 39, the CPU is unable to 
determine the correct sign because the mathematical methods cannot be used to extract the 
direction, either + a or - a, from the DOP data. As discussed above with respect to Fig. 9, 
compensating apparatus 700 determines the correct sign by measuring the DOPsign after the 

15 delay module. In other embodiments, referring to Fig. 10, the CPU of a compensating 
apparatus 800 determines the sign without monitoring a DOPsign value. Compensating 
apparatus 800 includes a fixed PMD module 810, such as a birefringent crystal or delay 
module shown in figure 3, which adds a fixed, known PMD time delay, Qknown, to the 
variable, unknown PMD time delay, ^unknown, produced within the optical fiber. By adding 

20 the fixed, known delay to the optical signal, CPU can identify the leading and trailing PSP by 
determining the correct sign in equation 39. CPU can determine the correct sign provided 
that Qknown is greater than Qunknown- For example, given a relationship such as z = 
(100+x) 2 , unambiguous values of x can be calculated for any measured value of z provided 
that x is less than 100. The exact value of Qknown depends upon the magnitude of 

25 ^unknown ^ can be estimated by a system operator prior to operation or can be determined 
by the CPU. For example, the CPU can multiply a previously calculated PMD time delay, 
^calc , by a scaling factor to produce a value of the fixed, known time delay that satisfies the 
relationship Qknown > ^calo The CPU, in turn, sends an electronic signal to fixed PMD 
module 810 to adjust the fixed delay to produce the necessary delay. 
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Referring to Fig. 11, a system 900 for calibrating a polarization transformer 940 
includes a transmitter 905, a polarization modulator 910, a polarimeter 920, a CPU 930, and 
a calibrating module 945. As described above, during operation, CPU 930 sends control 
signals to polarization transformer 940 to transform an optical signal such that a delay 

5 module (not shown) reduces PMD by retarding the leading PSR The objective of system 900 
is to determine, for a given set of electric control signals, what polarization state input into 
polarization transformer 940 will be transformed to different polarization states defined by 
the calibrating module. Calibrating module 945 includes a polarizing beam splitter (PBS) 
950 having an orientation p , and two photodetectors 960 and 970. As a transformed SOP, 

10 s , from polarization transformer 940 propagates into calibrating module 945, PBS 950 
redirects different portions of s onto detectors 960 and 970 based on the polarization of s . 
Each detector produces separate currents that are proportional to different portions of the 
transformed SOP. The currents are given by the expressions 

/i^Cl + s- p) (4°) 

15 and 

h = V 2 {\-s- p) (41) 
where h and h are normalized such that + i 2 ) = 1. An analog device 946 receives U and h 
and sends an electric signal proportional to the difference of the two currents, U, to CPU 930. 
id is the projection s • p . Alternatively, CPU 930 can receive h and i 2 and calculate the 

20 difference i d . As polarization modulator 9 1 0 changes the SOP of the optical signal 
transmitted by transmitter 905, CPU 930 determines the SOP of the optical signal by 
analyzing the electric signals from polarimeter 920, and records i d , and the magnitude of the 
control signals being sent to polarization transformer 940 for that SOP. By collecting a set of 
SOP measurements, Sj, and storing them into a matrix S. CPU can construct an 

25 overdetermined matrix given by 

S P = i d (42) 
where each Sj is a 3-dimensional row of S, p is a 3-dimensional vector, and i d is the vector 
of difference currents. CPU can estimate p by computing the pseudo-inverse of S. To 
calibrate polarization transformer 940, CPU computes an estimated p for a plurality of 

30 electric control signals and SOPs. CPU stores each estimated p , SOP, and the values of the 
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electric control signals in a lookup table. The lookup table provides specific electric control 
signals that will align the SOP with the axis system of the PBS. When compensating PMD, 
the calibrating module is replaced with the delay module such that the relative orientation 
between the module and PBS is known. After determining the location of the PSPs, the 

5 magnitude of the time delay, and which PSP is leading, the CPU uses the lookup table and 
the known relative orientation between the module and the PBS to determine which electric 
control signals are needed to transform the PSPs such that the delay module retards the 
leading PSP. In other embodiments, referring to Fig. 12, a calibrating module 945 and a 
delay module, such as described above with respect to figure 3, can be integrated into a 

1 0 single unit 95 1 so that the CPU can create a lookup table specifically for the polarization 
transformer and the delay module being used in a compensating apparatus. 

In other embodiments, the polarization transformer can include polarization 
transforming devices, e.g., electrooptic, acoustooptic, or stress induce bifringence, which can 
transform the PSPs to linearly polarized PSPs aligned with the x- and y-axes of the delay 

1 5 controller. The delay module also can include any optical device capable of delaying of one 
polarization state relative to another polarization state. Although shown aboveas separate 
units, the polarization transformer and delay module can be integrated into a single unit that 
functions to transform the optical signal and delay the leading PSP to reduce PMD in the 
optical signal. An integrated polarization transformer and delay module is described, for 

20 example, by C. Glingeneret. al. in Optical Fiber Communications Conference PD29, 1999. 
A measure of performance of the PMD compensating apparatuses can be determined by 
monitoring the jitter in the rising edge of the PMD compensated optical signal. A histogram 
of the rising edge includes a variance that is proportional to residual PMD. 

It should be understood that the foregoing description is intended to illustrate and not 

25 limit the scope of the invention, which is defined by the following claims. Other aspects, 
advantages, and modifications are within the scope of the following claims. 
What is claimed is: 
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